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Abstract 

The Higgs boson pair production process at hadron collider provides an opportu- 
nity for performing a study of the trilinear Higgs boson self-coupling. In this paper, 
we analyze the pair production of a neutral Higgs boson via both gluon-gluon and 
b-b fusions in the littlest Higgs (LH) model at the CERN LHC. We find that in some 
parameter space the relative corrections of the total cross section to the SM prediction 
may reach a value of 24% when x (— Afv'/v 2 ) — 0.95 at the LHC. We conclude that 
if the parameter x has a value above 0.8, the relative corrections contributed by the 
LH model reach values beyond 8% and can be observed at the LHC. 
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1 Introduction 

The standard model(SM) ^Hj theory has been proved by all existing precise experimental 
data with its theoretical predictions beyond one-loop level being coincident with exper- 
imental observations. But in the SM the Higgs boson mass suffers from an instability 
under radiative corrections. This "hierarchy problem" motivates much of current research 
works about new physics beyond the SM. Among the extended models beyond the SM, 
the little Higgs model offers a very promising solution to the hierarchy problem in which 
the Higgs boson is naturally light as a result of nonlinearly realized symmetry [3]-(H]. The 
first successful model, which cancels all relevant quadratic divergences based on the pseudo 
Goldstone idea, was constructed by Arkani-Hamed, Cohen, and Georgi [3]. Then more 
models were constructed, such as SU(5)/SO(5) 0, SU(Q)/SP(6) [7j, and the minimal 
moose SU(3) 2 /SU(3) [5j and general moose SU{3) n /SU{3) k 0. The most economical 
model of them is the littlest Higgs (LH) model, which is based on an SU(5)/SO(5) non- 
linear sigma model [5]. It consists of a SU(5) global symmetry, which is spontaneously 
broken down to SO(5) by a vacuum condensate /. In the LH model, a set of new heavy 
gauge bosons (Ah, Zh, Wh) and a new heavy- vector- like quark (T) are introduced which 
just cancel the quadratic divergence induced by SM gauge boson loops and the top quark 
loop, respectively. 

One of the most important task of present and future experiments is to search for 
Higgs boson and investigate its properties. Studying the properties of the Higgs potential 
will reveal details of mass-generation mechanism in spontaneously broken gauge theories, 
which can be obtained through measuring the Higgs boson self- interactions. Multiple Higgs 
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boson production processes at hadron colliders provide the way to probe the Higgs boson 
self- interactions. Many works have been contributed to studies of Higgs-pair production 
at the hadron collider in some traditional models ^Hj-EH- The possibility of measuring 
the Higgs boson self-coupling at the LHC has been investigated by U. Baur, T. Plehn and 
D. Rainwater [Tl|. They found that it should be possible at the CERN Large Hadron 
Collider (LHC) with design luminosity to establish that the SM Higgs boson has nonzero 
self- coupling and that X/Xsm can be restricted to a range of — 3.7 at 95% confidence 
level if its mass is between 150 to 200 GeV. Recently, the LH model contribution to the 
Higgs decay width was investigated in Refs. ^S] and ^H]- Dit> et al, discussed also the 
LH model contribution to the process pp — * H°H° + X in Ref. 15 . There they did 
not consider mixing and interference effects between the SM particles and the new heavy 
states, and thus they got negligible results of the order of (v//) 4 and concluded that the 
contribution from the LH model to the pair production to the Higgs bosons seems to be 
unobservable at the LHC ^3]. If the interference and mixing effects are included in the 
analysis, the contribution is at the order (v/f) 2 and does change obviously the results as 
compared to the SM prediction in some parameter space |16j . 

In this paper we investigate the effect of the LH model on neutral Higgs boson pair 
production via both gluon and bottom fusions at the CERN Large Hadron Collider (i.e., 
gg — * H°H° and bb — ► H°H° ) at the complete lowest order, considering mixing and 
interference effects between the SM particles and the new heavy states. The paper is 
organized as follows. In Sec. 2, we briefly go through the LH model theory. In Sec. 3, we 
present an analytical evaluation. The numerical results, discussion and a short summary 



3 



are given in Sec. 4. Finally we present the relevant Feynman rules in the Appendix. 



2 littlest Higgs model 

The littlest Higgs model is based on an SU(5)/SO(5) nonlinear sigma model. At the 
scale A s ~ 47r/, the vacuum expectation value (VEV) associated with the spontaneous 
symmetry breaking proportional to the scale / is parameterized by the 5x5 symmetry 
matrix [Hj ^Zl 



The VEV breaks the SU(5) global symmetry into its subgroup 50(5) and breaks the 
local gauge symmetry [577(2) (g) ?7(1)] 2 into its diagonal subgroup SU(2)l (g> U(l)y at the 
same time, which is identified as the SM electroweak gauge group. The scalar fields are 
parameterized by 



where H(x) = Tr a (x)X a is the Goldstone boson matrix. X a are the broken generators of 
SU(5) which obey the relation 




(1) 



S(a?) = e 



,ai(*)// s in(x) r // 



(2) 



X a S -S A ai =0. 



(3) 



Therefore, the Goldstone boson matrix n(a;) can be expressed as 

/ ht/y/2 <fi \ 





where 
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are a doublet and triplet under the unbroken SU (2) l®U(1)y SM gauge group, respectively. 

The leading-order dimension-2 term for the scalar fields in the littlest Higgs model 
can be written as 



is the covariant derivative for gauge group [SU(2) <g> U(l)] 2 = [SU(2)i U(l)i] 
[SU{2) 2 ® U(l) 2 ]. It is defined as 



where W w - = ^Li w &Qj and B J = B w Y 3 are the SU ( 2 )i and U ( 1 )j g au g e fields ' 
respectively. The generators of two ST/(2)'s (Qj) and two J7(l)'s generators (Yj) are as 
follows 



Y 1 = diag{-3, - 3, 2, 2, 2}/10, Y 2 = diag{-2, - 2, - 2, 3, 3}/10, (8) 



where (j a (a = 1,2,3) are the Pauli matrices. As we expect, the breaking of the gauge 
symmetry \SU{2) x U{\)] 2 into its diagonal subgroup SU(2)l x U(1)y gives rise to heavy 
gauge bosons W and B', and the remaining unbroken subgroup SU (2)lxU (l)y introduces 
the massless gauge bosons W and B. 

As we know, in the littlest Higgs model there is no Higgs potential at the tree level. 
Instead, the Higgs potential is generated at one loop and higher orders due to the interac- 
tions with gauge bosons and fermions. Up to operators of dimension 4, the Higgs potential 
(Coleman- Weinberg potential) can be expressed as HZ] 



(6) 
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V = A^ 2 / 2 Tr(0t0) + i\ h ^ h f(h^h T - h*<t>tf) - tfbtf + \ h 4hh*) 2 
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+XhWhhtf<f>tf + A^/i/JTr^V) + X^ 2 [Tr(0V)l 2 + A^Tr^Vf^), (9) 



where A^2, A/^/j, A^4, \h(f><f>hi ^h 2 4> 2 -> ^4> 2 <t> 2 , an d A^4 are the coefficients of the original Higgs 
potential. The coefficients which are concerned in our calculation have the expressions as 

m 

\ _M^ 2 _xM (j ? _M^ 2 4M/ 

By minimizing the Coleman- Weinberg potential, we obtain the vacuum expectation values 
(h°) = v/y2, (i<jP) = v' of the Higgs boson doublet h and triplet 4>, which give rise to 
electroweak symmetry breaking (EWSB). The Coleman- Weinberg potential provides the 
trilinear and quartic Higgs self-couplings. After EWSB, the gauge sector gets an additional 
mass and mixing term due to the VEVs of h and eft. By diagonalizing the quadratic term 
of the gauge sector, we may get the mass eigenstates Al, Zl, Wl, Ah, Zh, and Wh and 
their masses. 

To avoid large quadratic divergence in the Higgs boson mass due to the heavy top 
quark Yukawa interaction, we introduce a pair of new fermions t and t' |S] an d a set of 
new interactions. The scalar couplings to the top quark can be taken from the following 
Lagrangian [H] [T7] : 

C Y = ^\ife l]k e xyXl ^ j ^ k yu'i + X 2 ftt /C + H.C. (11) 

where x = (&3> *3) *)> e ijk and e xy are antisymmetric tensors where i, j, k run through 
1, 2, 3 and x, y run through 4, 5, and Ai and A2 are the new model parameters. By 
expanding the above Lagrangian, we get the physical states of the top quark t and a new 



heavy-vector-like quark T. The masses of the two physical states are 
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v\ 1 X 2 J u 

mt = 7Wm\ + T 2 



1 fv' 1 A? f, \\ 



(12) 



m T = f^\j + Xl[l + 0(v 2 /f 2 )], (13) 

respectively. Since the top quark mass is already obtained in the SM, we can then get the 
parameter relation from Eq. (fT2|) as deduced in Ref. JZJ 

1 1 v 2 

aTaTw* 2 ' <14> 

3 Calculation 

At hadron collider, the Higgs boson pair can be produced through two mechanisms. One is 
loop-induced production via gluon fusion; the other is from bb annihilation. The Feynman 
diagrams contributing to the subprocess gg — > H°H°, which are involved in the framework 
of the LH model, are depicted in Fig. 1. The diagrams created by exchanging two external 
gluon lines and final Higgs boson lines are not shown there. In Fig. 1, the first three 
diagrams(Fig. 1(1)-1(3)) [excluding the Fig. 1(1) with a heavy- vector-like quark (T) 
loop] are just the same as those in the framework of the SM, while the remaining four 
figures (Fig. 1(4)-1(7)) are the extra diagrams beyond the SM. All the Feynman diagrams 
can be classified into three types. The first type is named as s-channel diagrams with 
the exchange of a virtual neutral Higgs boson H° or heavy triplet Higgs boson <J>° which 
couples to a pair of gluons via a triangle quark loop [shown in Fig. 1(1), 1(4)]. The second 
type is called box diagrams [shown in Fig. 1(2), 1(3), 1(5), 1(6)] and the third is the 



quartic interaction type where the neutral Higgs bosons are produced by means of quartic 
interactions [shown in Fig. 1(7)]. All the relevant Feynman rules can be found in Ref. 
|17j and the Appendix in this paper. In the loop diagram calculation of this subprocess, 
we adopt a dimensional regularization scheme. The Feynman diagrams for the subprocess 
bb -> H°H° are depicted in Fig. 2. In this work we adopted the Feynman-'t Hooft gauge. 
In our calculation, we denote the two subprocesses as 

g(c u Pi)+g(c 2 , p 2 ) - H°(k 1 ) + H°(k 2 ) , (15) 
b(c 3 , pi) + 6(04, P2) -> H°{k 1 ) + H°(k 2 ) . (16) 

where c±, c 2 , C3, and C4 are the color indexes of the two initial particles (gg and 66). pi, p 2 
and ki, k 2 are the incoming and outgoing four-momenta of the initial and final particles, 
respectively. In the subprocess gg — ► H°H°, we define # as the scattering angle between 
one of the gluons and one of the final H° bosons, and in the subprocess bb — ► H°H°, 6 
represents the scattering angle between the b quark and one of the final Higgs bosons. In 
the center-of-mass system (c.m.s.) the four-momenta of the initial and final particles can 
be expressed as 1 



Pi 

P2 

h 



4 0, 0, f ), 



{—, —p H osw.0, 0, —f3 H ocos9), 



V S y S y S 

k 2 = (—, -— P H os'me, 0, - — f3 H o cosO), 



(17) 



In our calculation we set the mass of the b quark to be zero except in the Yukawa coupling. 



8 



where (3 H o = yl — Am? H0 /s is the velocity of final neutral Higgs bosons in the c.m.s. and 

S = (Pi +P2) 2 - 

The amplitude for the subprocess gg — > H°H° can be expressed as 
M(gg -»• H°H°) = M {t) + M {b) + M {q) 

^ClC2 

+he a ^k laP2(3 + he^ PlaPw ) , (18) 

where .M^, M^ b \ and .M^ represent the amplitudes of triangle, box, and quartic dia- 
grams, respectively. 

For the subprocess bb — > H°H°, the amplitude can be expressed as 

M(bb^H°H°) = v( P 2)(gi+92h)Sc 3 cMPi)- (19) 

In the above two equations /j (i = 1, 5) and gj (j = 1, 2) are the form factors of the two 
subprocesses, respectively. Since the explicit expressions of these form factors are lengthy, 
we do not list them in this paper. 2 

Then the total cross sections for these two subprocesses can be written as 

a(s, gg - H°H°) = J* d ?£ \ M(gg - H°H°) | 2 , (20) 

1 f i+ 

a(s, bb -► #°F°) = ^-^ y, dt J] | M(bb^ H°H°) | 2 , (21) 

respectively, where the bar over the summation recalls averaging over initial spins and 

colors, i* 1 = {m 2 H0 — 1/2 ± s/3 H o/2), and due to the identical final two Higgs bosons, the 

2 The Mathematica program codes of all the form factors for gg — > H°H° and bb — > H°H° are obtainable 
by sending an email to moonbank@mail.ustc.edu.cn 



right-hand sides of Eqs. Q2U|) and l|21jl have been multiplied by an additional factor of 1/2 
separately. The color and spin average factors for the subprocess gg — > H°H° are 1/64 
and 1/4 and for the subprocess bb — > H°H° are 1/9 and 1/4, respectively. The total cross 
section for the neutral Higgs boson pair production through gluon (or bottom) fusion in 
proton-proton collisions can be obtained by performing the following integration 



In our calculation we adopt the CTEQ6 parton distribution function ^S] and take the 
factorization scale \i to be 2m^o in the subprocess gg — > H°H°, while /i to be m H o/2 in 
the calculation of the subprocess bb — ► H°H° |19| . The numerical calculation is carried 
out for the LHC at an energy of 14 TeV. 

4 Numerical results and discussions 

In the numerical evolution we take the input parameters as m\y = 80.423 GeV, mz = 
91.1876 GeV, m u = 4.5 MeV, m d = 8.5 MeV, m s = 150 MeV, m c = 1.25 GeV, m t = 
174.3 GeV, and a(mz) = 1/128. We use a simple one- loop formula to express the running 
strong coupling constant a s : 




(22) 



where yfs and \/s denote the pp and gg (or bb) c.m.s. energies, respectively, and dLij/dr 



is the luminosity of colliding partons, which is defined as 




(23) 



a s (m z ) 



(24) 



1 + 



33-2n/ 

67T 



10 



where Oi s { m z) = 0.118 and nj is the number of active flavors at scale \i |2U| . 

In the numerical calculation, we use the next-to-leading order formula to evaluate the 
running mass of bottom quark fn b (Q) 21 . 



m b (Q) = U 5 (Q,m b )m b (m b ), Q < m t , 

m b (Q) = U e (Q,m t )U 5 (m t ,rn b )fn b (m b ), Q > m t , (25) 

where m b = m b (m b ) = 4.25 GeV and the energy scale Q is taken to be 2m H o in our 
calculation. The evolution factor Uf(f = 5,6) is 

\a s (Qi) J I 4vr 



32-2/' 
rf _ 8982 - 504/ + 40 f 

J ~ 3(33 -2/)2 ' [Zb) 

In the LH model, the calculation of the Higgs boson pair production at hadron colliders 
involves four additional free parameters. One is the parameter /, which is the symmetry 
breaking scale parameter at TeV order. The direct and indirect effects of the LH model 
provided by the present experimental measurements have placed a constraint / > 3.5 TeV, 
although it depends on the model assumption about the U(l)'s [12] |23| , The second 
one is the mass of the Higgs boson, rriffo. The third parameter is the coefficient A2, which 
is the coupling constant of the new heavy- vector-like quark T. Because there is a relation 
between the top quark and heavy- vector-like quark coupling parameters Ai and A2, as 
shown in Eq. (|14|). we can use A1/A2 to parametrize the mass of the new heavy vector 
quark [see Eq. (|13|) ]. The last one is v', and we define x = Afv'/v 2 to parametrize this 
vacuum expectation value of the scalar triplet field eft. The masses of neutral scalar boson 
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M$o can be given as ^5] UH 

2m\ f 2m 2 H0 f 2 
*° v 2 [l-(Av'f/v 2 ) 2 } v 2 (l-x 2 ) ' 1 ' 

The above equation about the mass of requires a constraint of < x < 1 (i.e., 4u' / jv 2 < 



1), which shows the relation between the scale / and the vacuum expectation values of 
the Higgs field doublet and triplet (v,v r ). Based on the limitation of current electroweak 
experimental data ^7] [22] [23], in our calculation we choose / = 3.5 TeV unless otherwise 
stated. 

Since our numerical calculation shows that the contributions of subprocess bb — ► H®H° 
to the parent process pp — ► H°H° + X at the LHC is less than 1% of the contributions 
from subprocess gg — > H°H°, we present only plots of the subprocess — ► H°H° for 
representation. In fact, our analysis also demonstrates that the contributions of the ad- 
ditional diagrams involving the new heavy-vector-like quark T or neutral heavy triplet 
Higgs boson $° for the subprocess gg — > H°H° are very small. The deviations of the cross 
section from the SM are mainly aroused by the contributions of the diagrams which exist 
also in the SM, but the interactions between Higgs bosons and quarks in the LH model 
are different with the corresponding ones in the SM. 

The relative effect of the LH on the cross section [5 = {(?lh — &sm) / &sm] for subprocess 
gg — > H°H° as functions of the c.m.s. energy (V§) of incoming gluons for the cases of x = 
(corresponding to M$o = 3.06 TeV) and x = l/\/2 (corresponding to M$o = 4.32 TeV) 
are depicted in Fig.3(a), with the parameters taken as / = 3.5 TeV and m H o = 150 GeV. 
In both figures, the solid line, dashed line and dotted line correspond to A1/A2 = 1/2, 
A1/A2 = 1, and A1/A2 = 2, respectively. From Eq. (|12 |l -l|14 |) we can get that those curves 
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correspond to tut = 6.28 TeV, 5 TeV, and 6.28 TeV, separately. All the curves in Fig. 
3(a) have the common line structure which decreases rapidly in the vicinity of the Higgs 
boson pair production threshold with an increment of V§, and then increases steadily after 
arriving at its minimal value. The curves also obviously show that with an increase of the 
value of A1/A2, the effect of the LH model is getting stronger. 

In order to clarify the line shape in Fig. 3(a) more clearly, we present the cross sections 
in the LH model (aLH), the SM (osm), and the difference between them (&lh — &sm) 
for the subprocess gg —> H°H° as functions of V§ with m H o = 150 GeV, / = 3.5 TeV, 
A1/A2 = 1, and x = in Fig. 3(b) and Table 1. The solid line, dashed line and dotted line 
in Fig. 3(b) correspond to &sm, and &lh — &SM, respectively. We can see from Fig. 
3(b) that in the VI region from 310 GeV to 500 GeV the curve for &sm rises up steeply, 
while the &lh — &sm changes gently, which can be also read from Table 1. Fig. 3(b) and 
Table. 1 show that the polelike behavior of the c.m.s. energy around \fl ~ 400-500 GeV 
in Fig. 3(a), comes from the fact that the cross section(<T,sAf ) rises up steeply when V§ 
is just beyond the threshold energy and decreases gently after &sm reaches its maximal 
value \/I ~ 500 GeV, while the variation of &lh — &sm is relatively slow in our plotted 
energy range. 

The dependence of the ratio o^hI^sm of the subprocess gg — ► H°H° on the parameter 
x with m#o = 150 GeV and \fl = 800 GeV, is depicted in Fig. 4(a). The solid line, dashed 
line and dotted line correspond to A1/A2 = 1/2, A1/A2 = 1, and A1/A2 = 2, respectively. 
Fig. 4(a) shows that the effect of the LH model always enhances the cross section of 
gg — ► H°H° in our chosen parameter space. The three curves also demonstrate that the 
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U sm 




f J j H ~ CSM 

i-J 1 A. t_7 1 VI 


WT,H ~ & SM ) /<5"<?M 


310 


0.0052071 


0.0060828 


0.0008756 


0.168162 


330 


0.0268346 


0.0295622 


0.0027275 


0.101643 


350 


0.0884609 


0.0947308 


0.0062699 


0.0708781 


370 


0.2052696 


0.2169084 


0.0116388 


0.0567002 


400 


0.3425389 


0.3595995 


0.0170606 


0.0498064 


420 


0.4051675 


0.4245253 


0.0193578 


0.0477772 


440 


0.4476802 


0.4685783 


0.0208981 


0.0466808 


460 


0.4738619 


0.4957345 


0.0218726 


0.0461582 


490 


0.4903020 


0.5128817 


0.0225796 


0.0460525 


500 


0.4912129 


0.5138794 


0.0226666 


0.0461441 


520 


0.4880278 


0.5107043 


0.0226765 


0.0464655 


540 


0.4797802 


0.5022973 


0.0225171 


0.0469321 


590 


0.4462748 


0.4679308 


0.0216559 


0.048526 


700 


0.3565471 


0.3754663 


0.0189191 


0.0530621 


900 


0.2334130 


0.2478193 


0.0144063 


0.06172 


1000 


0.1935886 


0.2063001 


0.0127115 


0.0656626 


2000 


0.0598705 


0.0659044 


0.0060339 


0.100783 



Table 1: The cross sections of gg — > H°H° in the LH model and SM, the difference between them, 
and the relative correction with m H o = 150 GeV, / = 3.5 TeV, Ai/A 2 = 1, and x — l/y/2 



effect of the LH model is not sensitive to x in the range of x < 0.8, but their values 
increase rapidly and can be larger than 1.1 when x > 0.85 for all three curves. To explain 
the result that the correction blows up as we take the x — > 1 limit shown in Fig. 4(a), we 
decompose the cross sections of subprocess gg — > H°H° in the LH model into three parts: 



&LH = &box + &tri + Pint, (28) 



where we denote at, ox , atri, and &i n t as the contribution parts from the box diagrams [including 
quartic diagrams, shown in Fig. 1(2), 1(3), 1(5), 1(7)], the triangle diagrams [shown in Fig. 
1(1), 1(4)] and the interference between the box and triangle diagrams respectively. In Fig. 
4(b), we show the contributions of these three parts as functions of x with the conditions 
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of y/§ = 800 GeV, m H o = 150 GeV, / = 3.5 TeV and \ 1 /X 2 = 1. The solid line, dotted 
line, dash dotted line, and dashed line correspond to box diagrams, triangle diagrams, 
interference contributions, and total cross section, respectively. Our calculation demon- 
strates that the main contribution to the cross section of subprocess gg — > H°H° comes 
from the box diagrams, and the interference contribution is mainly from the Feynman 
diagrams involving the top quark. Let us review the couplings of H°-H°-H° and H°-i-t 
in the LH model. They can be expressed as 

. / 3m 2 33m 2 v x 2 \ , . 

9hoh°h° : -(^f-- ^ x _ x ^ (29) 



9H°tt : — 
v 



xv 2 x 2 v 2 2 v 2 v 2 \{ ( \\ 



2/ 2 4/ 2 3/2 /2 A 2 + A 2 ^ Al 2 + A 2 2 



• (30) 



We can see that the contribution of the triangle Feynman diagrams is s-channel suppressed 
and is relative small, due to the heavy 3> and a/I = 800 GeV >> m H o. The figure 
shows that the contribution from the interference between the box diagrams and triangle 
diagrams is negative and blows up quickly when we take the x — ► 1 limit. We also find 
that the dependence of the product of gH°H°H° an d 9H°tt on the parameter x behaves with 
the same rapid increment when x is close to 1. Therefore, we can conclude that the quick 
enhancement behavior of the total cross section in the vicinity where x — > 1 arises mainly 
from contributions of the interference between the triangle diagrams involving the gnoH H 
vertex and the box diagrams involving the g h 0f t coupling. As we know, perturbativity alone 
should put some bounds on the range that x should be allowed to vary. But with these 
limitations some couplings in the original Higgs potential must behave badly. Since we 
started with a relatively will-behaved Higgs potential, it is clear that x — > 1 cannot be a 
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very well-defined limit and should not be considered as a physical limit. In other words, the 
original Higgs potential is an effective potential result of integrating out the heavy states. 
Although it cannot give definitive limitations, the size of all Higgs potential parameters 
should be roughly order 1 theoretically. We depict the relations between x (or qhhh) and 
the absolute values of original Higgs potential parameters A^2, A/^, A^, and \h<f>4,h [see 
Eq. (|10|) ] in Fig. 4(c) in the conditions of / = 3.5 TeV and m#o = 150 GeV. In the figure, 
the solid line, dashed line, dotted line, and dash-dotted line correspond to | A^2 1 , IA/^1, 
\\h±\, and \\h(j><t>h\-> respectively. The absolute value of the coupling gH0H°H° corresponding 
to x value is scaled on the upper axis. From the figure, we can see that the values of all 
four Higgs potential parameters, which are concerned in our calculation, are far beyond 
order 1 when x > 0.95 and |A/^0/i| is less than 10 _1 when x < 0.25. Therefore, the above 
x value range should be excluded in our consideration, since these absolute values of the 
Higgs potential parameters cannot satisfy the request of the well-behaving effective Higgs 
potential. Then we can approximately put bounds o the varying range of the parameter 
x as 0.25 < x < 0.95 on the condition of / = 3.5 TeV and m^o = 150 GeV. 

The ratio of total cross sections ulh/^sm^ as functions of x at the LHC (y/s = 14 TeV) 
with / = 3.5 TeV and A1/A2 = 2 is depicted in Fig. 5. The total cross section involves the 
contributions of both subprocesses gg — ► H°H° and bb — ► H°H°. In the figure, the solid 
line, dashed line, and dotted line correspond to m H o = 115 GeV, 150 GeV, and 180 GeV, 
respectively. We can see again that the deviation of the total cross section in the LH 
model from the corresponding SM value is not sensitive to x when x < 0.8, but increases 
quickly when x > 0.9, and the relative deviation, defined as S = {(Jlh — &sm)/&sm-, can be 
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beyond 12% when x reach 0.9 at the LHC. We can see from this figure that the correction 
effect of the LH model is obviously related to the Higgs boson mass. When the Higgs 
boson mass varies in the range of 115 — 180 GeV, the curves show that the heavier the 
Higgs boson mass is, the stronger the correction effect becomes. 

In Fig. 6 we plot the ratio between the total cross sections in the LH model and the 
SM, as the functions of the mass of the neutral Higgs boson m^o with the conditions of 
/ = 3.5 TeV and A1/A2 = 2 at the LHC. In this figure, the solid line, dashed line, and 
dotted line correspond to the cases of x = 0.25, l/y/2 and 0.9 respectively. The relative 
correction of the LH model to the SM cross section [5 = (am — 0sm)/(tsm\, can reach 
14.8% when x = 0.9 at the LHC. 

The total cross sections asM an d am of process pp — > H°H° as functions of the mass 
of the neutral Higgs boson m H o is depicted in Fig. 7 at the LHC, with / = 3.5 TeV, 
x = 0.9, and A1/A2 = 2 in the LH model. The solid line and dashed line correspond to 
the cases of the SM and LH model, respectively. They show that the absolute correction 
induced by the LH model decreases with an increase of the Higgs boson mass. 

In Fig. 8, we plot the total cross sections am of process pp — ► H°H° as functions 
of the new heavy vector-like quark mass my with m^o = 150 GeV and x = l/\/2 in the 
LH model at the LHC. In this figure the solid line is for A1/A2 = 1/2, dashed line for 
A1/A2 = 1, and dotted line for A1/A2 = 2. According to Eq. (|13|) for a fixed value of 
A1/A2, the mass of heavy- vector- like quark T is only related to the symmetry breaking 
scale parameter /. For example, when A1/A2 = 2 and tut varies from 1.5 TeV to 4.0 TeV, 
we have that the parameter / changes from 0.85 TeV to 2.26 TeV. We can see from this 
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figure that in the chosen parameter space the cross sections are enhanced when either m-r 
is relatively small or A1/A2 is relatively larger. 

In conclusion, we investigated the effect of the LH model on the pair production process 
of neutral Higgs bosons via both gluon and bottom fusions at the LHC. The numerical 
analysis shows that with the possible parameters, the relative cross section correction to 
the SM prediction may reach a value of 24% at the LHC when x = 0.95. We conclude 
that when the parameter x has a value above 0.8, the relative correction contributed by 
the LH model reached a value beyond 8% and could be observable at the LHC. 
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Foundation of China and a grant from the University of Science and Technology of China. 

Appendix 

The interactive Lagrangian of the scalar field S and the up-type quarks of the first two 
generations take the same form as in Eq. (Jll[). except that there is no need for the extra 
t. The interactive Lagrangian of the S and the down-type quarks can be expressed as 



where the isospin index i = 1,2. All the Yukawa couplings between the Higgs boson and 
quarks (except the Yukawa coupling which involves the top quark) can be obtained from 
the corresponding Lagrangian directly. We present the expressions of the H°-H°-u-u and 
H°-H°-d-d couplings in the equations 



£y = -XdftijktxyXi^jx^tyd 



+ H.C. 



(31) 



Amdv' 



4m u v' 



(32) 



9H°H°dd 
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9H°H°uu 
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The expressions of other couplings concerned with in this work can be found in Ref. |17j . 

In the LH model, the trilinear interaction of Higgs bosons gH°H°H° S e ^ s a correction 
to the SM at the order of v 2 / f 2 , and an additional H°H°& coupling is generated. They 
are given by the Lagrangians 

-Luhh = H°H°H° (- 2fV ' Xh * h + vA h , - 12 V '" Xhi + 6 ^ XhHh 

\ V V V 



-C HH4> = H°HV + 14V2 " f f Xh * h - 6V2 v'X hi + ^p^j , (33) 



From the above Lagrangians we obtain the Feynman rule for H°H°H° and H°H°&° 
couplings as [T7j 

3m| 66M| t/ 2 



9h°h°h° '■ -i To 

l y fly 



. /56V2M|>' 3 2V2Ml v' 29y/2 M| ?/\ 
: -* I -4 ~2 I • ( 34 ) 
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Figure Captions 

Fig. 1 The one-loop Feynman diagrams of the subprocess gg — > H°H° in the LH 
model: (1), (4) s-channel diagrams. (2), (3), (5), (6) box diagrams. (7) quartic diagrams. 
The notations u and d represent up-type and down-type quarks, respectively. The dia- 
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grams created by exchanging two external gluon lines and final Higgs boson lines are not 
shown. 

Fig. 2 The lowest-order Feynman diagrams of the subprocess bb — ► H°H° in the LH 
model. 

Fig.3(a) The relative effect of the LH on the cross section [S = (&lh — &sm) /&sm\ for 
subprocess gg — > H°H° as functions of y/i on the conditions of x = \/y/2J = 3.5 TeV, 
ruffo = 150 GeV. The solid line, dashed line, and dotted line correspond to A1/A2 = 1/2, 
A1/A2 = 1, and A1/A2 = 2, respectively, (b) The cross sections of gg — > H°H° subprocess 
in the SM, the LH model, and their difference cflh—vsm as functions of y/1 with x = 1/ y/2, 
f = 3.5 TeV, m H o = 150 GeV, and Ai/A 2 = 1. 

Fig. 4(a) The dependence of the ratio olhI&sm for the subprocess gg — ► H°H° on 
the parameter x with m#o = 150 GeV and y/i = 800 GeV. The solid line, dashed line and 
dotted line correspond to A1/A2 = 1/2, A1/A2 = 1, and A1/A2 = 2, respectively. (b)The 
total cross section and contributions of the box diagrams (including quartic diagrams), 
the triangle diagrams, interference between the box and triangle diagrams (&lh = &box + 
&tri + &int) for the subprocess gg — > H°H°, as functions of x with the conditions of 
y/I = 800 GeV, m H o = 150 GeV, / = 3.5 TeV, and Ai/A 2 = 1. The solid line, dotted 
line, dash dotted line, and dashed line correspond to box diagrams, triangle diagrams, 
interference contributions, and total cross section respectively. (c)The absolute value of 
the original Higgs potential parameters A^2, Xhcf>h, A/j4 and \h<f><f>h as functions of x with the 
conditions of m H o = 150 GeV and / = 3.5 TeV. The absolute value of coupling gH°H°H° 
corresponding to the x value is scaled on the upper axis. 
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Fig. 5 The ratio of the total cross sections in the LH model and SM (Jlh/&sm for 
the process pp — ► H°H° at the LHC, as functions of x (= Afv'/v 2 ), with / = 3.5 TeV and 
Ai/A 2 = 2. 

Fig. 6 The ratio of the total cross sections olhI^sm for the process pp — > H°H° at 
LHC, as functions of the Higgs boson mass m H o with the conditions of / = 3.5 TeV and 
A1/A2 = 2 in the LH model. 

Fig. 7 The total cross sections asM and clh for the process pp — > H°H° at the LHC, 
as functions of the Higgs boson mass m H o with / = 3.5 TeV, x = 0.9 and A1/A2 = 2. 

Fig. 8 The total cross section glh as functions of the mass of the new heavy- vector-like 
quark itlt with m H o = 150 GeV and x = l/\/2 in the LH model at the LHC. 
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